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Abstract 

The paper is studying the dynamics of twisted vortex bundles, which were detected in exper- 
'^■» 
O 
/*> iniental investigations of superfluid turbulence in superfluid '^He-B. The analysis shows that a 

^O linear torsion oscillation of a vortex bundle is a particular case of the slow vortex mode related 

(N 

with the inertial wave, which was already investigated in the past in connection with observation 

(D of the Tkachenko waves in superfluid He and the experiments on the slow vortex relaxation in 

Q superfluid ^He-B. The paper addresses also a twisted vortex bundle terminating at a lateral wall 

Cd of a container starting from the elementary case when the bundle reduces to a single vortex. The 

s 

I theory considers the laminar regime of the vortex-bundle evolution and investigates the Glaberson- 

^H Johnson-Ostermeier instability of the laminar regime, which is a precursor for the transition to 

O 
I ^1 the turbulent regime at strong twist of the bundle. The propagation and the rotation velocities of 

fSJ the vortex front (the segment of the vortex bundle diverging to the wall) can be found from the 

> 

OO equations of balance for the linear and the angular momenta, and the energy. It is demonstrated 

^^ that the vortex front can move with finite velocity even in the absence of mutual friction (the T = 

t^^ limit). The theory is compared with experimental results on vortex- front propagation in superfluid 

3He-B. 
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I. INTRODUCTION 

Theoretical and experimental investigations of various vortex structures in classical and 
superfluid hydrodynamics have a long history^ ^. In particular, helix vortices, which ap- 
pear in wakes of propellers and other spinning bodies, were intensively studied in classical 
hydro dynamicaSI. Recently, they became an object of investigation also in superfluid hydro- 
dynamics, however, not as individual single vortices, but as bundles of helix vortices, also 
called twisted vortex bundles. The interest to twisted vortex bundles arose in connection 
with experimental studies of the transient process of establishing of stable vorticity in ro- 
tating superfluid ^He-BP'. The twisted vortex bundle appears in an originally vortex-free 
rotating container with a superfluid when vortex lines being injected at container's bottom 
expand into the rest part of the container. The experimental and theoretical studies of the 
twisted vortex bundle addressed possible turbulence of the vortex front and the transition 
from the laminar to the turbulent regim#"^. Meanwhile a proper dynamical theory of the 
laminar regime of the propagating vortex front is still lacking. There were estimations of 
the connection between the dissipation rate and the velocity of the vortex front, but the 
front velocity was not derived from the vortex dynamics. However, without a satisfactory 
theory of the laminar regime it is impossible to have a reliable physical picture of the more 
complicated turbulent regime or of the transition between the laminar and the turbulent 
regime. The first step in this direction was undertaken in Ref. [HI but for the simplest case 
when the vortex bundle is not twisted and the force driving the vortex front is very weak. 
The goal of the present work is to suggest a more realistic and general theory. 

The paper starts from the overview of the uniformly twisted vortex bundle in Sec. |TTj This 
is a generalization of the previous analysis^ on the case when the vortex-line tension is impor- 
tant. The section presents the analysis of the Glaberson-Johnson-Ostermeier instability^^Hs] 



of the twisted vortex bundle (Sec. II C), which demonstrates that the laminar regime be- 
comes unstable at rather weak twist. The instability can be considered as a precursor of the 
transition to the turbulent regime. Section |III| analyzes linear oscillations of the vertically 
uniform vortex bundle on the basis of the general theory of slow vortex oscillations developed 
decades agcP'. Section IV presents the dynamical theory of a vortex bundle terminating at 



a lateral wall starting from the elementary case when the vortex bundle reduces to a single 
vortex. The balance equation for the linear and the angular momenta, and the energy are 



derived and then used for calculation of the vortex front velocities of propagation along and 
rotation around the container axis. A special attention is devoted to the limit of vanish- 
ing mutual friction (the T = limit) and to qualitative comparison with the experiment. 
Concluding discussion is presented in the last Sec. IVl 
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FIG. 1. Untwisted and twisted vortex bundle, (a) Untwisted vortex bundle. Any vortex line is in 
an axial plane, i.e., is not twisted. Such a structure corresponds to the equilibrium state in the 
coordinate frame rotating with the angular velocity il (Ref. [TT1) . (b) Twisted vortex bundle. The 
twisted stem of the bundle is at rest in the coordinate frame rotating with the angular velocity O. 
The vortex front (whorl) propagates with the velocity Vf along the axis z and rotates around the 
axis z with the relative velocity Ail = Qf — 0, with respect to the container, the normal fluid, and 
the bundle stem rotating together. 



II. UNIFORMLY TWISTED STATIONARY VORTEX BUNDLE 

A. Continuous vorticity 

We start from the purely continuum approach to superfluid vorticity neglecting quantum 
elastic effects (line tension of vortex lines and Tkachenko shear rigidity), which are propor- 
tional to the circulation quantum n. So vorticity inside an untwisted vortex bundle of radius 
i?o in a container of radius R [Fig. [l]^a)] is continuous and uniform, the azimuthal superfluid 
velocity field being 

J Qr r <Ro, 
\^-^ R> r > Ro 

In a uniformly twisted vortex bundle [Fig. [lib)] any vortex line forms a helix of pitch 
271 /Q moving along a cylinder coaxial with the z axis, where Q = wq^/wo^t = u^/uzT is the 
wavenumber of the twist (called later twist), uq^ and uqz are the azimuthal and the axial 
components of the vorticity cjq = V x 1?^ , i;^ is the superfluid velocity, and u^ and Uz are the 
azimuthal and the axial components of the displacements u of vortex lines from positions in 
the untwisted bundle described by Eq. ([I]). In an untwisted bundle the vorticity ljq has the 
only axial component Uq = ujqz = 2Q. In a twisted bundle azimuthal components Uq^ and u^ 
appear, but radial components uor and Ur still vanish. In the continuum approach the twist 
Q{r) can be an arbitrary function of the distance r from the axis. This is a consequence 
of ignoring Tkachenko shear rigidity of the crystalline vortex bundle: shear deformation of 
the bundle does not cost any energy. On the other hand, choosing a r-independent twist 
Q one may expect that Tkachenko rigidity will not have an essential effect on the twisted 



bundle. Later it will be shown (Sec. IIIB) that this expectation is true only for a stationary 
twist, while considering the torsional oscillation of the bundle the role of Tkachenko rigidity 
is important. 

We are looking for the stationary twisted state in the coordinate frame rotating with the 
angular velocity Q. In the continuum theory the vortex line moves with the velocity vl 
coinciding with average superfluid velocity Vs (Helmholtz's theorem), and the vortex lines 
are at rest if the superfluid velocity in the rotating frame has no component normal to the 
vortex line. This yields the conditiorpl 

[vs^ir) - nr]uJozir) - Vszir)uJo^{r) = 0. (2) 



One can check by substitution that the following pattern satisfies this condition: 
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r dr (l + gV2)2' "-P dr {1 + Q^r^y ' 

where Vq is a constant determined by a natural condition that the total axial mass fiow 
27r J Vsz{'r')rdr must vanish in a closed container. In the vortex-free region R > r > Rq 
the axial velocity Vsz does not depend on r and ii R ^ Rq Vsz must be very small. This 
means that the absence of the axial mass flow requires that Vsz{Rq) = and correspondingly 
Vo = QQR^. In the opposite limit Ro ^ R the bundle occupies the whole cross-section of 
the container, and the axial mass flow vanishes if 

Q^R' 



Vo 



n 
Q 



_ln(l + g2/?2) 

Then the velocity field in the twisted bundle is 
Qr Q^R^ Q 



-1 



Vs<p 



2d2 



Q'R 



g2r2) ln(l + g2/^2) 



-1 



(4) 



(5) 



1 + g2r2 ln(l + g2i?2) ' g 

Knowing this field one can find the energy e and the angular momentum m^ per unit length 
of the bundle. For the bundle filling the whole container {R ^ Rq) they are given by 



2nps 



"^ v% + viz 
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(6) 



(7) 



Twisting of the bundle leads to the fiux of the angular moment along the vertical axis. Using 
the known hydrodynamic expression 



Uij = P6ij + PsVsiV. 



Jsj 



(8) 



for the momentum fiux tensor of the superfiuid component, the fiux of the angular momen 
tum along the axis z is given by 



Jrn = 2-nps I Vs^VszT dr 
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Q' 



.In2(l + g2i?2)(l + g2i?2) 



(9) 



The proper thermodynamic potential for a liquid in a container rotating with constant 
angular velocity ^2 is the Gibbs potential g = e — Qrriz- An important property of the 



continuous- vorticity approach is that the Gibbs potential of the uniform vortex bundle is 
always negative independently from how strong the twist Q is. But in the presence of 
quantum-vorticity effects, namely, line tension, the condition g = Q can be satisfied as 
shown in the next subsection. This important for dynamics of the vortex-front propagation 



(Sec. IV) 



B. Quantum-vorticity effect: vortex-line tension 



Now we generalize the analysis of the twisted vortex bundle by taking into account 
vortex-line tension. The widely accepted method of dealing with the line-tension effects is 
the Hall-Vinen-Bekarevich-Khalatnikov (HVBK) theoryHMl This theory uses the local- 
induction approximation, when the contribution of line tension to the energy and the su- 
perfluid velocity depends on the circulation quantum and the curvature radius of the vortex 
lines. Helmholtz's theorem is still valid but now the vortex velocity vl coincides with the 
local superfluid velocity Vsu which differs from the average superfluid velocity Vs and in the 
local-induction approximation is given by 



v,i = v, + z/^V X s. 



(10) 



where s = loq/uq is the unit vector tangent to the vortex lines, Us = (k/Att) ln(r^/rc) is 
the line-tension parameter, and r^ and Vc are the intervortex distance and the core radius 
respectively. Then Eq. ([2]), which provides the condition that vortex lines are at the rest in 
the rotating coordinate frame, transforms to 

z/,Q(2 + QV2) 



^s,ir) - fir + ^^"'^,",^,1, 



ujQzir) - 



r + 



Wo0(r) = 0. 



;i + QV2)3/2 _ 
The superfluid velocity and the vorticity fields satisfying this condition are [cf. Eq. (|3 
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1 + gV ' (l+QV2)3/2' 
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i^sQ' 



^0()> 
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(12) 



:i + gV2)2 ' (l + Q2^2)5/2 ' -^<P (^ + Q2^2)2 ' (^ + Q2^2)5/2 

Further in this subsection only the case R^ Rq is considered. The axial mass flow vanishes 
if 
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(13) 



The density of the Gibbs potential g = e — Qniz includes the vortex line-tension term 
proportional to Ug-. 

"R rR 



g = 2ttPs 



^ vL + vt 



rdr + 2tt Pal's / i^o^ df — 2,1x0. ps 



v.^r dr. 



(14) 



Substituting the superfluid velocity field given by Eq. (12) and bearing in mind that the 



modulus of vorticity vector is Uq = Uqz a/1 + Q'^R'^ one obtains that 



n ^n p2^o nps{n + QvoY 
g = TvpsK — 
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(15) 



g2 \^^ (i + g2i^2)i/2^ 

Line tension contributes not only to the energy, but also to the momentum fiux tensor, 
which becomes^^'^ 



liij = P5ij + PsVsiVsj + Ps^s <^o5ij - 






(16) 



Using this expression one can find the axial fiux of the angular momentum along the axis z: 

' r dr. 



Jrt 



i-R pR 

2ttPs Uszt'^ dr = 2nps I \ v sci,v sz - - s- 

Jo Jo V ^0 



ly^ 



(17) 



After substitution of the values of the components of the velocity Vg and the vorticity ojq 
and following integration over the cross section of the cylindric container one obtains 



Jyy 



Q' 
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where 



Vs = Vs 



2 + Q^R^ 



-2 



(1^ 



(19) 



/I + g2i?2 

is the renormalized line-tension parameter taking into account the bundle twist Q. The 
expression for the fiux of the angular momentum along the axis z can be also derived from 
the thermodynamic definition Jm = ~de/dV z'^\.m^ = —dg/dVz^\^, where subscripts mean 
that derivatives with respect to Q = Vzf are taken at fixed rriz or ^2 respectively. But 
this method is complicated by the condition that the derivatives must be calculated at 
constant number of vortices, whereas derivation of the fiux of the angular momentum from 
the hydrodynamic momentum-fiux tensor is more straightforward. 



The line tension terms oc Ug make possible the condition g = 0, which provides the balance 
of forces on the vortex front separating the vortex bundle from the vortex-free region. In the 
limit QR — )■ the condition g = determines the angular velocity Q = Sz/^/i?^, at which the 
vortex bundle terminated at the lateral wall is in the equilibrium with the rotating container. 
This velocity was found in Ref. [Hi In the opposite limit QR — ?■ oo the condition is satisfied 
aX Q = 2vsQ/ R- So the strong twist Q can facilitate the balance of forces on the vortex front 



and make possible the stationary motion of the vortex front without friction (see Sec. IV). 



C. Stability of twisted vortex bundle 

An important feature of a twisted vortex bundle is a mass flow along vortex lines. Long 
ago Glaberson et al^^^ showed that this may course instability with respect to excitation 
of vortex array waves. Let us discuss the threshold for the Glaberson- Johnson-Ostermeier 
instability starting from the case of a uniform superfluid rotating with the angular velocity 
VL. 

The general spectrum of the linear plane waves propagating in a rotating superfluid 

^Y pipz+ik-r—iujt : J14I 

-' = (2^ + -sP') (2^^ + -sP' + 1^) , (20) 

where p is the wavenumber along the z axis, k is the wave vector in the xy plane, and 
ct = \J kVl/Sti is the Tkachenko-wave velocity determined by the shear rigidity of the 
vortex lattice. Further in this subsection we shall neglect the Tkachenko contribution since 
the calculation has shown that its effect is not essential. The Landau critical velocity for 
the superflow along the vortex line for instability with respect to excitation of a wave with 
wave vector (fc,p) is 



v.{k) = ^ = ^(2fi + u,p-) (^ + u^ . (21) 

The minimum of this velocity at small p <^ a/2^/V^ and large k ^ y^2n/us is the critical 
velocity for the Glaberson- Johnson-Ostermeier instability: 



Vg 



^2VtVs. (22) 



So the instability starts not from pure Kelvin waves {p ^ k) but from waves with p <^ 
k propagating normally to the rotation axis. Still it would be interesting to discuss the 
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instability for pure Kelvin waves k = when Vc = {2Q + VsP'^)/p- Analyzing the stability 
of a single vortex in a rotating coordinate frame Ostermeier and Glabersonl^ obtained 
Vc= (fi + VsV^^IV- The factor 2 difference in the limit p — )■ is due to long-range interaction 
between vortices. In contrast to the case of a single vortex, the HVBK theory of the vortex 
array considers a uniform displacement of the whole vortex array. So any vortex feels the 
velocity induced not only by its own displacement but also by displacements of other vortices. 
This is valid until the wavenumber p does not exceed the inverse intervortex distance l/r^, ~ 
^/Vl/k. If pr^ ^ 1 the velocity field induced by an oscillating vortex is cut off at the distance 
of the order the wavelength 27r/p, and the gap of the Kelvin-wave spectrum becomes Vt as 
in the single-vortex case. The difference between the Kelvin-wave spectra for pr^ ^ 1 and 
pTy <ti 1 was noticed by Raja Gopal^^ long ago. 

Switching to instability in a twisted vortex bundle one should replace the uniform vorticity 
2VL in the expression for vq by the space-dependent vorticity modulus uJoir), i.e., vq = a/wqi^. 



This approach is justified since a;o(r) varies slowly at the scale ^Ji7j2VL at which instability 
sets on. The Glaberson-Johnson-Ostermeier critical velocity vg should be compared with 
the superfluid velocity Vi = {v^ — Qr)s^+VzSz along the vortex lines in the rotating coordinate 
frame. For a large number vortices when QR^ ^ Us one can use calculations for continuous 



vorticity (Sec. II A). Then the twisted vortex bundle is stable only for rather weak twist 



satisfying the inequality 



Q^-n^ < j^. (23) 



III. TORSIONAL OSCILLATIONS OF THE VORTEX BUNDLE 

A. Phenomenological approach 

One should expect that weak stationary uniform twisting is related with some soft (Gold- 
stone) torsion mode related with axial symmetry of the vortex bundle. Let us start from 
the phenomenological approach to torsional oscillation: there are weak oscillations of the 
angular momentum density accompanied with the oscillating flux of the angular momen- 
tum. Restricting ourselves with a single degree of freedom related to a pair of canonically 
conjugated variables "2; component of the moment m^ -the angle ip of rotation around the 



axis z" the phenomenological Hamiltonian is 

where I is the moment of inertia and A is torsion stiffness. The Hamihon equations are 

dip 5% mz{z) 
dt Sniz I 

drriz ST-L „ , ,„ , ,^^. 

V. (AV.y^) . (25) 



dt 5ip 

For a plane wave ex: e^^^"*"^* propagating along the axis z these equations yield the dispersion 
relation 

u^' = vy = ^p\ (26) 



where v^ = yA/I is the velocity of the torsion mode. Equation (25) confirms that the 
quantity Jm = —AVz^p can be considered as the angular momentum flux along the z axis. 
The effective Hamiltonian "H can be derived from the Gibbs potential of the liquid with 
the vortex bundle. Keeping in mind that Q = V^v^ the torsion stiffness is determined by 
expansion of g with respect to Q. In the continuum model, Sec. II A this gives A = 
'kPsQ'^Rq/3 for the thin bundle i?o < -R and A = irp^Q'^R^/U for the thick bundle R ^ Rq. 
As for the inertial term it follows from the expansion of the Gibbs potential g{^') with 
respect to i7' — i7 = d(p/dt. This yields the moment of inertia I = 7tPsRq/2 independently 



of the ratio R/Rq. Then according to Eq. (26) the oscillation dispersion relation is 



for Rq <^ R and 



u' = In'Ry (27) 



u^ = In^Ry (28) 
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for Rq ^ R. 

B. Analysis based on the linear vortex dynamics 

The simple analysis given above relied on the assumption that even in the excited state the 
bundle rotates as an ideal solid body with the velocity field given by Eq. ([I|. The assumption 
is not self-evident, and moreover, not always valid as demonstrated below. A more accurate 
approach must start from the vortex dynamics equations with proper boundary conditions. 

10 



The linear hydrodynainic equations for a superfluid with vorticity were presented in Ref. 
[TH For slow motion described by the plane wave oc e«'«''+«p^-«'^* these equations are: 



2 



^^ icok 

- lUVsW - 2ilVL± —Vs\\ = 



-ioovs± + 2nvL\i = 0' (29) 



C'-prC 



vl± = -iujux ^ Vs±, (30) 

where u and vl = —iuou are the displacement and the velocity of the vortex lines from their 
equilibrium positions in the xy plane, and subscripts || and ± denote vector components in 
the xy plane parallel and normal to the in-plane wave vector k. The equations correspond 



to the dispersion relation given by Eq. (20) with the line-tension terms oc Vg neglected: 

^^' = 40^7^ + 4^', (31) 

The slow mode under consideration is a combination of the inertial wave (the first term 
oc VL^) and the Tkachenko wave (the second term oc c\). This combined mode (called the 
mixed mode in Ref. |T3) was crucial for interpretation of experiments on Tkachenko waves 
in ^He and on the slow mode in ^He-B (see Refs. [T^ and [T7|) . 

The torsional oscillation mode is a plane wave along the z axis and an axisymmetric 
cylindric wave in the xy plane. In the configurational space in the cylindric system of 
coordinate the equations of motion for the axisymmetric slow mode (no dependence on the 
azimuthal angle 0) are 



dz"^ dz"^ V dr"^ r dr r^ 

du 



- 2f] Vif + I ^^^ + -^^ - ^ ) = 0' 

+ 2ilUr = 0, 



dt 

dt "^ 2Vl \ dr'^ r dr r"^ J ' 
where radial (subscript r) and azimuthal (subscript 0) components correspond to the lon- 



gitudinal (subscript ||) and transverse (subscript ±) components in Eqs. (29) and Eqs. (30) 



respectively. The general solution of these equation finite at the axis of the container is 

u^ = u^Ji{kr)e'^'-'^\ Vsr='^^Jiikr)e'^'-'^\ (33) 
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where Ji{kr) are the Bessel functions of the first kind. The axial velocity 

^^^^ moJo(A;r)e*P^-*-* (34) 



is determined from the incompressibility condition 



^ + ^ + ^ = 0. (35) 

dz dr r 

The next step is to formulate proper boundary conditions. The first one is imposed on the 

liquid velocity and follows from the analysis of the vortex-free region, where the velocity is 

divergence- and curl-free. The general axisymmetric velocity field satisfying these conditions 

is 

v,r{r) = [Ah{pr-) - BK,{pr)]e'P\ Vs^ = 0, t;,,(r) = i[AIo{pr) + B Koipr)]e'^' . (36) 

Here Io{pr) and KQ{pr) are the modified Bessel functions of the first and the second kind. 
The constant A and B are determined by the boundary condition a.t r = R: 

Vsr{R) = [AhipR) - BK,{pR)]e'P' = 0, (37) 

and by the continuity condition on the boundary of the vortex bundle r = Rq: 

vURo) = [AhipRo) - BKi{pRo)]e'P\ 

v^Ro) = ^[AIo{pRo) + BK,{pRo)Y^\ (38) 



where Vsr{Ro) and Vsz{Ro) are the superfiuid velocity components determined by Eqs. (33) 



and (34). 



At pR ^ 1 this imposes the following condition on the velocity at the boundary of the 
vortex bundle: 

ivURo) - vURo f^^l'^'^ = 0. (39) 

For a thin bundle R ^ Rq this reduces to the condition Vsz{Ro) = 0, while if the vortex 
bundle fills the whole container {R — Rq <^ Rq), the radial velocity component must vanish: 
VsriRo) = 0. This reduces to the condition 27r J^ '^ Vsz^dr = of zero axial mass flow used 
for the stationary twisted state. 

The second boundary condition is imposed on the shear components of the vortex-lattice 
stress tensor, which must vanish at the boundary of the vortex bundled 



CfltriRo) = —PsCx 



dr 
12 



= 0. (40) 

r=Ro 



First let us look for a solution of the boundary problem neglecting the Tkachenko shear 



rigidity {ct — )■ 0). Then the boundary condition (40) is not relevant. Let us focus on the 



case of a thin vortex bundle Rq <^ R. Other cases can be treated similarly. li Rq '^ R the 
boundary condition Vsz{Ro) = requires that JoikRo) = 0. The slow mode corresponds to 
the lowest nonzero root of the Bessel function: k = 2A05/Ro- This yields the frequency of 
the slow torsion mode: 

k^ 



u 



2 -^ 4^2^ 



2 0.692n^Rlp^. 



(41) 



Though this frequency only slightly differs from the frequency in Eq. (27) (numerical factor 
0.692 vs. factor 0.667) obtained from the phenomenological approach the origin of the 
difference is worth of discussion. The solution of dynamical equations does not support the 
assumption of the phenomenological approach that the vortex bundle oscillates as an ideally 
solid body. Indeed, Eq. (33) shows that the azimuthal velocity Vs^j, = —iuu^ is proportional 



to the Bessel function Ji{kr) where k ~ I/-R0 aiid is different from the solid-body rotation 
when Vs(j) oc r. On the other hand, the question arises, whether ignoring of Tkachenko 



rigidity is justified in this case. If p is small in the dispersion relation (31) the Tkachenko 



term c^k"^ ~ kQ/R^ definitely exceeds the other term related with the classical inertial wave 
and therefore cannot be neglected. 

Thus we must take into account the Tkachenko rigidity and look for possible in-plane 



wavenumbers k satisfying the dispersion relation (31). The latter is a bi-quadratic equation 
for k so two values of k"^ are possible: 



ki 



1 



CO 



-p 




+p' - 



4^2 



.pz 



(42) 



Since uj ~ QRop is proportional to p [see Eqs. (27) and (28)] but much larger than ctP this 
expression at small p reduces to 



kl 



\c? 2VL 

2 c^ Ct 



(43) 



It is important that |A;-|-| is much larger than p but still much smaller than the inverse radius 
1/ Rq. The latter condition will allow to expand Bessel functions (see below). The general 
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solution of the boundary problem is a superposition of two modes with k^ and A;_: 



Ji{k+r) Ji{k_r) 



VLt 



du. 



Ua 



k\ +p2 
k+Jo{k+r 



dr 



u+Ji{k+r)+u^Ji{k_r)], Vsz = i2np -+ , 2 , n 

I k^ + p"^ 

dv> 
u+k+J2{k+r) - U-k-J2{k_r), u^ = -— ^ 



u. 



k_ +p2 
k_Jo{k_r) 



k_ +p2 



r 



dv.. 



dr 



2ilipu^. (44) 



The constants u+ and U- are determined by the boundary conditions (39) and (40). In 



the case i?o ^ -R the determinant of these two linear equations vanishes at the condition 

(kl + p^)k+J2ik+Ro) {kl +p^)k.J2ik.Ro) 



k+Jo{k+Ro) 
4(1 + klRl/Q) - kt{l + klRl/Q) 



k.Joik.Ro) 



u 



2^2 






0. 



(45) 



Here the expression (43) for k^ and the condition k±Ro <^ 1 were used. Equation (45) yields 



exactly the same dispersion relation (27) as the phenomenological approach. It is interesting 



that this dispersion relation does not contain the circulation quantum, which determines the 
Tkachenko rigidity. Despite it the rigidity is important since it provides a solid body motion 
of the vortex bundle. 

The analysis of the present section ignored mutual friction focusing on the T = limit. 
Meanwhile, whatever weak mutual friction is, at very long wavelength of the torsion (slow) 
mode its spectrum transforms from sound-like to diffusive. For the clamped regime when 
the normal component co-rotates with the container as a single solid body (the case relevant 
for superfiuid ^He) this was demonstrated in Sec. VII. F of Ref. [TH The diffusive slow mode 
in ^He-B was investigated experimentally and theoretically in Ref. [T71 In connection with 
the torsional oscillation of the vortex bundle the diffusive slow mode was discussed by Eltsov 
et a/.H. 

IV. THE BUNDLE TERMINATING AT THE WALL: PROPAGATION OF THE 
VORTEX FRONT 

A. Single vortex line terminating at the lateral wall 

Before considering a vortex bundle terminating at the lateral wall (the main topic of this 
section, see Fig. [I]) it is instructive to analyze a much simpler problem of a single vortex 
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terminating at the lateral wall. This case allows an analytical solution, which helps to 
understand what may be going on in the vortex bundle. 

While without mutual friction the vortex line moves with the local velocity of the su- 



perfluid Vgi (see Sec. II B), in the presence of mutual friction the velocity of the vortex line 



depends on the mutual friction parameters a and a': 

vl = Vsi + a[s X {vn - Vsi)] - a[s x[sx {vn - Vsl)]], (46) 

where s is the unit vector tangent to the vortex line. Studying the effect of mutual friction 
one must know the external force on the whole liquid (per unit length of a vortex line): 

Psk[s X {vl - v,i)\ = psK{-d[s x[s X {vl - v„)]] + d'[s X {vl - v„)]}. (47) 

Here the other friction parameters d and d' were introduced, which were used by Bevan 

et al. ^^ and are more convenient for the present analysis. They are connected with a and 

a' by the relations 

d I 1 — d' , , 

"=rf2 + (l_rf/)2' ^~" =rf2 + (l_^/)2- (48) 

The vortex line goes along the container axis (the z axis) deviating from it continuously 
starting from some height and eventually terminates at the lateral wall. The end segment 
of the vortex line connecting the axis and the lateral wall can be also considered as an 
elementary case of the vortex front. Stationary propagation of this vortex front along the z 
axis supposes that the vortex moves as a solid body with constant vertical velocity f /, at the 
same time rotating around the z axis with angular velocity Vtf different from the container 
angular velocity VL. The normal liquid corotates with the container. This means that the 
relative velocity v^ — Vn inside the front is fully defined and has only two components: the 
z component Vf and the azimuthal component [Vlj — VL)r. 

For a single vortex in the local-induction approximation the superfluid velocity is deter- 
mined by the line-tension force: 

[s X Vsi] = -VsN. (49) 

Here N = ds/dl is the curvature vector with its magnitude equal to the inverse curvature 
radius. Two functions z{r) and 0(r) determine the shape of the vortex line in cylindric 



coordinates r,(l), z. Using Eq. (49 ) and the expression for the curvature vector in the cylindric 



coordinate frame the vector equation (47) leads to two equations for the axial and the 
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azimuthal components: 

Sr -{i^fr + u. 



Sr < -Vf + 



dz/dr 



dr [1 + r^{d(j)/dry + [dz/dr)^]^/^ 
Us d r'^{d(f)/dr) 



'l—s1)dvf+(srd'—SzSff,d){Vlf—Vl)r, (50) 



r dr [1 + r'^{d(j)/drY + {dz/dr) 



— I = (1 - sl)d{Qf - Q)r - [srd' + s,s^d)vf. 

(51) 



The third radial component of the force balance equation (47) is not an independent equation 



being a consequence of Eqs. (50) and (51) 



Integrating Eqs. (50) and (51) over the whole vortex line (keeping in mind that the line 



length element is dl = dr/sj. = ^yl + r'^^dcjy/drY + (dz/drydr) and taking into account the 
boundary conditions dz/dr = oo at r = and dz/dr = dip/dr = at r = i? one obtains 



VL 






d^VfR + {d' - d^^){VLf - Vt) 



R^ R^ 



R^ 
2 

i?2 



(52) 



(53) 



where three new mutual friction parameters related to the dissipative force oc d were intro- 
duced: 



d. 



d 
R 



R 



I -St 



dr, dd 



3d 



R 



l-sl 



-r dr, d. 



2d 



R 



°zS(f> 



r dr. 



(54) 



Equation (52) is the balance of axial forces on the vortex, while Eq. (53) is the balance of 



moments around the z axis. This becomes evident if one rewrites them as 



e — QfUiz = —psi^ 



dzVfR + [d! - dz^){Qf - Q)- 



VfTUz = -psK. 



p3 p2 

d^{nf -n)— -{d' + dz^)vf — 



(55) 



(56) 



where e = p^KZ/s is the energy and mz is the z angular momentum per unit length of straight 
vortex line far below the termination point respectively. For the vortex coaxial with the 
container and at the distance ri from the container axis 

R^-rl 



rriz = PsK,- 



(57) 



For the axial vortex considered here ri = 0. The left-hand side of Eq. (55) is a mechanical 



force on the vortex front, which is balanced by the axial friction force on the right-hand 



side. Equation (56) demonstrates the balance between the moment transferred to the liquid 
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because of the vortex propagation (the left-hand side) and the mutual friction torque (the 
right-hand side). All terms in the balance equations except for the dissipative terms oc d 
do not depend on the shape of the vortex. In particular, the terms oc Vf and Vtf originate 



from the Magnus- force term psn[s x Vl\ in Eq. (47). Both the terms contain the angular 
momentum of the single vortex state per unit length. 

The equations of the linear and and the angular momenta balance lead to the balance of 
the total energy, which determines the energy dissipation: 

dE 
~dt 



-psK 



d;,viR + d^iVtf - 9) 



Vf(e — VLfniz) + (fi/ — VL)vfm 
i?3 



dz^iytf — Vt)vfR 



(58) 



Solving Eqs. (52) and (53) one obtains 



2d^R{Q - Qo) 



Ad,d^ + 3[(1 - d'y 
3(1 - d' - d,^)in 



fin) 



(59) 



Adj^ + ?,[{i-d'Y-di^y 

where VLq = 2i>s/R^ is the critical angular velocity at which the Gibbs potentials of the single 
vortex state below the vortex front and of the vortex-free state above the front are equal. 
One cannot use these expressions directly for determination of the propagation and the 
rotation velocities vj and Qf because the mutual friction parameters d^, d^, and dz^ depend 



on the vortex shape determined from the differential equations (50) and (51). However 
assuming that the friction force is weak, either because d and d' are small or because the 
difference fi — fio is small one can determine the values dz, d^, and dz^, using the shape of the 
vortex line in the state of equilibrium solid-body rotation together with the container and 
the normal liquid [Eq. (13) in Ref. |TT]. The equilibrium assumes that Vf = 0, Q = Qf = Qq, 



and the vortex line lies in the axial plane: 

1 ry/2W 



s^ = 0, 



a/1 + (dz/dr) 



R^ 



VT 



R' 



i?2 



(60) 



Then 



d,. 



d^ = 3d 



d I p 

'0 

pdp 



2 - p'^dp 



2v^ 



-d = 0.609d, 



3(^2-1)^= 1.24ci, dz 



0. 



(61) 
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Using these values in Eq. (59) one obtains 

2A8dR{n - Qo) 






3.02^2 + 3(1 -dO^' 
3(l-c;')(^-^o) 
3.02^2 + 3(1 -rf')2^ 



the dissipation rate being 



dE __ d4,{n - no^R^ _ _ i.2Ad{n - noyn^ 

III ~ '^'"^Ad^d^ + 3(1 - d'Y ~ "^^''3.02^2 + 3(1 -rf')'' 
The general expressions for the dissipation rate and the vortex front velocity Vf [Eqs. 



(58) and (59)] point out that these quantities are determined not by the mutual friction 
parameter a only, contrary to previous estimations in the literature^, and in general the 
another parameter a' also affects the results. However, in the limit of the weak friction force 
one may neglect a tiny difference between the factors 3 and 3.02 in the denominators of Eqs. 



(62) and (63). Then according to Eq. (48) the dissipation rate and the vortex front depend 



only on a. In particular, in this limit the vortex front velocity is 

w/ ^0.83ai?((]-fio)- (64) 

Recently the single-vortex front dynamics was numerically simulated on the basis of the 
Bio-Savart law^. The results are in a qualitative agreement with those obtained here using 
the local-induction approximation. In particular, expansion of the expression (11) of Ref. [19] 
for the front velocity (noted as vlz there) with respect to i7 — r2o gives Vf ~ 0.74lQ;-R(fi — i^o), 



which does not differ essentially from Eq. (64). A more quantitative comparison requires 
more numerical data for lower angular velocities fi ~ fig while available dat#^ focus on high 
angular velocities fi :» ^q. 

Another useful approximation is to assume that the shape of the end vortex segment in 
the front is close to a horizontal straight line between the axis (r = 0) and the wall {r = R). 
Then s^ = s^j, = and s,. = 1, which corresponds to c/^^ = and d^ = d^ = d. In this case 
the final expression describing the vortex motion are 

2dR{n - fio) 






4rf2 + 3(l-rf')^' 

3(l-rf')(^-^o) 
4rf2 + 3(i_rf')2 ' 



dE _ djn - npfR^ 

dt ~ P''' M^ + 3{1 - d'r ^^^^ 
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The difference with the exact solution for the weak-friction case is not so pronounced, and 
we shall exploit this approximation for a more complicated case of the vortex bundle. 

B. Propagation of the vortex front: the energy and moment balance equations 

Let us consider now a bundle terminating at the lateral wall. The part of the bundle 
diverging to the wall is a vortex front separating the vortex-filled and the vortex-free parts 
of the container (Fig. [I]). Motion of the front along the container axis is a transient process 
of vorticity penetration into a container in the spin-up experiments^. Sometimes one can 
find stationary states of the bundle when the front does not move along the z axis and the 
bundle and the front (whorl) rotate as a solid body with the angular velocity determined 
from the thermodynamic analysis^. Here we address the case of a moving front. 

The front motion leads to change of the energy and the angular momentum and is accom- 
panied by mutual friction with the normal component moving rigidly with walls, or surface 
pinning and friction of vortex ends at rough surfaces of walls. We look for the state with 
the vortex front moving with constant velocity f j, while the angular momentum necessary 
for this motion being supplied via the vortex pinning at the bottom and transferred to the 
front by the flux of the angular momentum J^,. The flux is related with the twist Q = Vzf 
of the vortex-bundle stem. Let us start from observation that a uniformly twisted bundle 
stem in the wake of the front moving with the velocity Vf is possible only if the front rotates 
with respect to the bundle stem with the relative angular velocity 

An = nf-n = VfV^^. (66) 

Only this purely kinematic relation provides that the moving front leaves behind it a uni- 
formly twisted bundle in the absence of vortex reconnections. 

The rotation and the vertical motion of the front with respect to the solid-body rotating 
container with the normal fluid lead to friction. The balance of the angular momentum 
around the axis z during the front motion is 

m^Vf + Jm= Tfr. (67) 

The equation tells that the angular momentum brought by the angular-momentum flux 
■Jm = —de/dVzflfn^ ~ ~^5'/^^2V^|q ^^ compensated by the growth of the total angular 
momentum (the term nizVf) due to front propagation and by the friction torque Tfr- 
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The balance of the hnear moinentuin along the axis 2;, as we shall demonstrate below, is 



e — JljTTlj. — W z^Jn 



Ffr. 



(68) 



The left-hand side of the equation is a driving mechanical force on the front balanced by 
the friction force Ff,.. Note that elongation of the bundle due to front propagation does not 
lead to variation of the linear momentum, because the latter is exactly zero above and below 
the front. So the mechanical force does not contain a term proportional to Vf similar to the 



term rrizVf in the balance equation (67) for the angular momentum. 



It is important also to consider the balance of the energy: 



(e — Qmz)vf = — (f^/ — Q)Tfr — VfFfr. 



(69) 



The left-hand side is the rate of the energy variation due to front propagation whereas the 
right-hand side is the dissipation rate. The Gibbs potential density g = e — Qniz, which 



appears in the left-hand side of Eq. (69), has also the dimensionality of a force and may be 



called effective force. It differs from the mechanical force in the left-hand side of Eq. (68) 



because the vortex front not only moves along the container axis, but also rotates. For 
an untwisted bundle {V z'-P = 0) studied in Ref. |T1J the two forces coincide. The balance 



equations (67) and (68) for the angular and linear momenta of the vortex bundle differ from 



the similar balance equations (55) and (56) for the single vortex by the presence of terms 



containing the angular momentum flux Jm- 



The balance equations (67) and (69) for the angular momentum and the energy seem 



evident. On the other hand, Eq. (68) of the linear momentum balance needs a justification. 



The simplest argument in its favor follows from the fact the three balance equations (67)- 



(69) are not independent, and it is easy to check that Eq. (68) follows from Eqs. (67) and 



(69). If one believes the two, one must accept the third one. It is possible also to derive 



Eqs. (68) and (67) directly from the equations of hydrodynamics. 



All terms on the leftt-hand side of the equations, i.e, those not connected with mutual 
friction, do not depend on the velocity and vorticity distribution inside the vortex front, but 
only on the distribution inside the bundle stem far below the front. Let us demonstrate this 
for terms originating from the Magnus-force term psK,[s x vl\ in Eq. (47). In Sec. IV A it was 



shown that the contributions of this term to the linear and the angular momentum balance 
of a single vortex are flfrriz and —Vfiriz respectively, where rUz is the z angular momentum 
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per unit length given by Eq. (57). For the vortex bundle one should find the summary effect 
of all vortices in the bundle. This requires integration of the vortex distribution over the 
bundle cross-section far below the vortex front. In particular, the Magnus force contribution 
to the linear momentum balance equation is 



Jo 2 Jq 

Here m^ is the angular momentum per unit length of the vortex bundle far below the vortex 
front, and the relation 

.„„(.) = 1*^ (70) 

r dr 

between the vortex density n^ and the azimuthal superfiuid velocity Vg^ was used. Similarly 
one can check that the Magnus force contribution to the angular momentum balance equation 

is —VflTlz. 

However, in order to derive the whole balance equations it is more convenient to start 
from the original hydrodynamic Euler equation for the superfiuid component. For the sake 
of simplicity we restrict the derivation with the continuous-vorticity limit since its general- 
ization on the case with line-tension effect and the mutual friction force is straightforward. 

Then the Euler equation is 

Ov 

-^ + {vs,V,)v, + V/i = 0. (71) 

The r and z components of this vector equation in the cylindric coordinates (r, 0, z) for the 
axisymmetric velocity field (no dependence on 0) are 

dVsr , dVsr , OVsr ^l^ , <9/i ^ ,_^. 

-j^ + ^sT^^+v^— ^ + ^— = 0, 72) 

ot or oz r or 

dvsz . dvsz dvsz . dfi , 

-^r+^sr^ + Vz^ + ^ = 0. (73) 

ot or oz oz 

For a stationary moving vortex front dvgr/dt = —Vfdvsr/dz and dvgz/dt = —Vfdvgz/dz. 

Choosing the chemical potential fi in the vortex-free region far above the vortex front as 



a reference point, i.e., assuming that /x = there, one can integrate Eq. (73) for the z 
component along the container axis r = 0. This yields the chemical potential at the container 
axis far below the vortex front {z = — oo): 

//(O, -oo) = J dz \Vf-^ - 2 ^7 ) = -VfVsziO, -oo) . (74) 

21 



As a next step one can calculate dependence of fi{r) on the distance r from the axis below 



the vortex front (further the argument z = — oo is omitted) integrating Eq. (72): 



IJ,{r) = VfVsz{0) - 



2 +X "■*''''>- 



(75) 



Now one can determine the total variation of the linear momentum of the liquid integrating 
the momentum flux component Ilzz = P + v1^ [Eq. (pi)] over the cross-section of the bundle 
stem bearing in mind that this is a momentum only of the superfluid component and the 
pressure is P = psli: 



27r 



R 



Il,,r dr = TiR^ 



VfVsziO) - 



vUO) 



+ TT / v'^^r dr + 2n 



ttR' 



VfVsziO) 



^L(o) 



R rR 

2 

sz 

Jo 

^ r^ dr 

2 ^ j„ , „d2 / „.2 "' 



f „ 



dri 



+ 271 I v^^r dr + -kR- / v 



L"'o ^1 

R 



r dr 



r 



TT / vl^rdr.ilQ) 
'0 Jo ' Jo 

Calculating integrals for the velocity field in the twisted bundle in the continuous-vorticity 

limit [Eqs. (tsl) and (111)] and using the kinematic relation (66) one can check that Eq. (76) 



yields exactly the left-hand side of the linear momentum balance equation (68). 



Next let us calculate the total friction force and the total friction torque integrating the 



mutual friction force per unit volume [cf. Eq. (47) 



ffr = dujQ {{vl - Vn) - s[s ■ {vl - «„)]} + d'uo[s X {vl - fn)]- 



(77) 



over the whole vortex liquid: 



/oo pn 

dz / ujq {(i [f/(l — si) — SzS^{Vlf — f2)r] + d'sr{^f — fi)r} r dr 
-oo Jo 

= 2npsK, / dz < Uod [f/(l — si) — SzSif,{Qf — Q)r'j — d' ^ {VLf — VL)r > r dr 

= dvfnip + d'{Qf — Q)mz, 

/oo pR 

dz / Uq {(i [{^f — il)r{l — s^) — SzS^Vf] — d'srVj^ r^ dr 
oo Jo 

= d{Qf — Q)mT — d'vfmz- (78) 



Here Uq is the absolute value of the vorticity vector ljq = 'V x Vg and 

}r)r dr, 



r-OO pR 

mp = 2TrpsK I dz a;o(l — ^1 — s 

-oo Jo 
"OO pR / 

rriT = 2npsK / dz Wq I 1 — s^ — 
-oo Jo \ 



2 ^-^^\r^dr 



Qr 



(79) 
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are the effective moments related with the dissipative parameter d. Only the vortex front 
region contributes to the bulk integrals, since there is no vorticity in the vortex-free region 
above the front, while below the front the relative velocity v^ — Vn has no component normal 
to the vortex lines and the integrand also vanishes. 

The dissipation rate depends only on the dissipative mutual friction parameter d: 

{Vlf - n)Tfr + VfFfr = d[vjmF + i^f - fi)^mT] = dyjimp + V^ip'^rriT). (80) 

The explicit expressions for the friction force and the friction torque allow to derive from 
the balance equations the relations determining the front velocity Vf and the twist Q = V^v^: 

\ — d' 1 — a' rriTQ j mp + Q'^rriT nizVL — e , , 

d a rriz "^ {nizVt — e)mz^ d(mp + rrixQ'^) 

As in the single-vortex case, the dissipative forces depend on the velocity field inside the 
front, and in order to find the moments rrip and mj- we approximate the vortex line in the 
front by straight line segments normal to the axis. This means that s^ = 1 and s^ = s^ = 0, 



and Wo = —dvcp/dz. After integration over z Eq. (79) reduces to 



-R pR 



pn pn 

mp = 2iT / Vs(j){r)r dr, mp = 2tc / Vs4,{r)r^ dr. (82) 

Jo Jo 

The further analysis restricts itself with the case of the bundle occupying the whole 
cross-section of the container {R ^ Rq). Let us neglect first the line-tension effect using the 



relations given in Sec. II A In the limit of strong dissipative friction parameter d compared 



to 1 — d' Vf and Q are expected to be small and the expansion in Q reduces Eq. (81) to 



1-a' 56 ^„ SnR 

> —QR, Vf -^ — -. 83) 

One may compare this result with experimental measurements of the twist Q and the 
vortex front velocity Vf at temperatures higher than about 0.5 Tc. Figures 15 and 17 in Ref. [9] 
show that Q grows and Vf drops with decreasing temperature. This qualitatively agrees with 



Eq. (83 ) since both d and {l — a')/a = {l — d')/d grow with decreasing temperature according 
to measurements of Bevan et al)^ for ^He-B at T > 0.6 Tc (see their Figs. 5, 6, 8, and 9). 
The results of the experimenlF and the numerical calculations of the HVBH equations^ at 
high temperatures were described by the expression Vf ~ aQR. It is worthwhile of noting 
that at q; :^ 1 — a' there is no difference between a and 1/d, so the quoted result differs 
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from Eq. (83) only by the factor 3/8. It maybe explained by inaccuracy of our assumption 
on the shape of the bundle inside the vortex front. 

From the position of the nowadays discussions of the T = limit it is interesting to 
consider the opposite case of vanishing mutual friction rf, d! — )■ 0. In this limit both the front 
velocity Vf and Q grow and 

I -a' 2^„ 3fi dVLR 

Thus with vanishing mutual friction the twist become extremely large, while the vortex 
front velocity vj vanishes in this limit. At the same time the relative angular velocity 
Vtf — Vt = Qvf remains finite and the front rotation is the most important source of energy 
dissipation. However, one should remember that friction will not fully disappear even in the 
T = limit: surface friction of vortex ends moving along a rough wall would restrict the 
velocity of the vortex front, as mutual friction does at T > 0. Even a more serious problem 
with this limit is the Glaberson-Johnson-Ostermeier instability of the laminar regime at 



strong twist demonstrated in Sec. II C 



According to Sec. II B the joint effect of line tension and twist drives the system close to 



the state where g = e — ilniz vanishes. If the system reaches this state Eq. (81 ) has a solution 
without mutual friction with the twist Q determined from the condition g = e — Qniz = 
and with the front velocity Vf = Jm/fn^. For large numbers of vortices the state is possible 



only for large QR, and using the expressions of Sec. II B in this limit Eq. (81) yields: 



The Glaberson-Johnson-Ostermeier instability at strong twist also puts possibility to reach 
this state in the laminar regime in question. However, this option can revive above the 
instability threshold, where strong oscillations can increase the energy e allowing to reach 
the condition e — firriz = at larger vortex number and weaker twists (see discussion in the 
next section). 

V. DISCUSSION AND CONCLUSIONS 

This work analyzed dynamics of twisted vortex bundles in rotating superfluids revealed 
in recent spin-up experiments on superfluid 3jjg_g|I19iiQj ^g ^^^e first step the linear dynamics 
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was investigated, which demonstrated that the torsion oscillation mode involving weak bun- 
dle twisting is a particular case of the linear vortex dynamics of the slow vortex mode. This 
mode was investigated in the past in connection with observation of the Tkachenko wave 
in superfluid ^He (Ref. [Ti| and the experiments on the slow vortex relaxation in superfluid 
^He-B (Ref. |l7j). The strongly twisted bundle was also investigated, but it was demonstrated 
that the Glaberson-Johnson-Ostermeier instability prevents reaching the strong-twist case 
in the laminar regime. 

The analysis addressed also a twisted vortex bundle terminating at a container lateral 
wall. The segment of the bundle diverging to the wall (vortex front, or whorl) is propagating 
along and is rotating around the container axis. The analysis starts from the case of a single 
vortex terminating at the wall, which allows an analytic solution for a weak force driving 
the vortex along the container axis. The three equations for the balance of the linear and 
angular momenta and for the energy were derived, which were further generalized on the 
case of the vortex bundle. It was supposed that the vortex front propagation occurs in the 
laminar regime, without vortex reconnections. This provides a kinematic relation between 
the linear velocity and the angular velocity of the front and allows to find the parameters 
of the vortex-front motion. In this equations only the dissipative components of the mutual 
friction force and torque require knowledge of the velocity and the vorticity distribution 
inside the vortex front, and their contributions to the balance equations were estimated 
using simple assumptions on the vortex line shape inside the front. 

The analysis of the vortex-front propagation was performed under the assumption that 
the vortex bundle below the front fills the whole container cross-section. But vortex-line 
tension always leads to the vortex-free region near the container lateral walP^nS g^f^ h^q 
equilibrium. So the analysis was done for the number of vortices exceeding the equilibrium 
vortex number. In the experiments on the vortex-front propagation one may expect that 
the number of vortices in the bundle is determined not by the condition of equilibrium but 
the conditions of creation of the vortex bundle at the container bottom, and this number 
could be more or less than the equilibrium value. Knowing this number one can generalize 
the present analysis including the effect of the vortex-free region. This effect was neglected 
for the sake of simplicity, but taking it into account would result in more accurate numerical 
coefficients in our final expressions. 

Another essential assumption of the analysis was strong pinning of the vortex bundle 
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at the bottom, which provided the constant angular-momentum flux towards the vortex 
front. Though pinning at the lower end of the bundle is quite possible, especially if the 
"bottom" is in fact an interface separating the B phase from the A phase with the periodical 
vortex structure, one may address also the case without pinning. Then coupling between 
the vortex bundle and the container can be provided only by mutual friction. This would 
slightly complicate the analysis since the angular-momentum flux proportional to the twist 
will be not uniform along the container axis. The numerical calculations by Eltsovl^, who 
neglected pinning at the bottom, showed that the twist varied behind the vortex front rather 
slowly, nevertheless. Then one may use the present analysis assuming that the bundle twist 
in the theory is the twist well behind the vortex front but still rather far from the bottom. 
The important question, which attracted great attention in the literature on superfiuid 
vortex dynamics nowadays, is: What is dynamics of the vortex front in the limit of vanishing 
mutual friction (the T = limit)? Physical and numerical experiments provided evidence 
that at low temperatures the laminar regime transforms to the turbulent regime. The 
Glaberson-Johnson-Ostermeier instability of the laminar regime at growing twist can be 
considered as a precursor of this transformation. The instability condition is determined by 
the twist Q whereas earlier they supposecF'^ that the transition to turbulence is governed 
by the ratio 1 — a'/a, which is the quality factor for Kelvin waves. However, the twist Q is 
connected with 1 — a'/a, and using this connection [Eq. ([83))] the laminar-regime stability 



condition given by the inequality (23) can be rewritten as 



The right-hand side of the inequality is rather small at a large number of vortices. So the 
Glaberson-Johnson-Ostermeier threshold predicts small values of 1 — a' /a at the transition, 
while in the experiments evidences for turbulence appear at 1 — a' /a of order unity or 
more^. Probably this means that there is an intermediate stage between the laminar and 
turbulent regime characterized by large vortex array fluctuations but still without essential 
number of reconnections. Previously the effect of reconnections was investigated numerically. 
According to Hosio et a/.™ reconnections, which were registered in numerical experiments, 
become essential only below 0.3Tc. Meanwhile, the laminar regime apparently becomes 
unstable at higher temperatures. 

The analysis of the intermediate regime is beyond the scope of the present work restricted 
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with the laminar regime. But one may speculate what results of the present analysis obtained 
for the laminar regime could be retained in the new regime. Since reconnections are still 



not essential the kinematic relation (66) remains in force. This validates our treatment of 
the balance equations, but with energy, forces and torques recalculated taking into account 
large fluctuations triggered by the Glaberson-Johnson-Ostermeier instability. It is natural 
to expect that fluctuations increase the energy of the bundle behind the front. This would 
facilitate reaching the condition e — Qm^ = 0, which allows the vortex front motion without 
friction at T — )■ at vortex numbers larger than those in the laminar regime. If it were 
realized it would explain observed saturation of the vortex-front velocity dependence at low 
temperatures. 
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